Linear fractional stable motion: a wavelet estimator of the ol 

parameter 



Antoine Ayache 
UMR CNRS 8524, Laboratoire Paul Painleve, Bat. M2 
Universite Lille 1 
59655 Villeneuve d'Ascq Cedex, France 
E-mail: Antoine . AyacheOmath . univ-lillel . f r 



Linear fractional stable motion, denoted by {^_f/.Q(i)}teE, is one of the most classical stable 
processes; it depends on two parameters H G (0, 1) and a G (0, 2). The parameter H characterizes 
the self-similarity property of while the parameter a governs the tail heaviness of 

its finite dimensional distributions; throughout our article we assume that the latter distributions 
are symmetric, that H > 1/a and that H is known. We show that, on the interval [0, 1], the 
asymptotic behaviour of the maximum, at a given scale j, of absolute values of the wavelet 
coefficients of {XH.a{t)}t^^, is of the same order as 2"-''^^"-'^/"); then we derive from this result a 
strongly consistent (i.e. almost surely convergent) statistical estimator for the parameter a. 

Key words: stable stochastic processes; statistical inference; wavelet coefficients; Holder regularity. 

1 Introduction and statement of the main results 

Let H and a be two parameters such that a € (1, 2) and 1/ a < H < 1. We denote by {XH,a{t)}t&. 
the symmetric a stable linear fractional stable motion (Ifsm for brevity) (see e.g. [51 [3]), defined, for 
all t € M, as, 
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Abstract 




(1.1) 



where (•) is a symmetric a-stable random measure and, for each z G M 



max{z, 0}. 



(1.2) 
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The parameter H characterizes the self-similarity property of Ifsm; namely, for all fixed positive real- 
number o, the processes {X{at)}t£R and {a^X(t)}tgK have the same finite dimensional distributions. 
The parameter a governs the tail heaviness of the latter distributions. The process {XH,ait)}teR 
has a modification with continuous nowhere differentiable sample paths; it is identified with this 
modification in all the sequel. 

The statistical problem of the estimation of H has already been studied in several articles: [6l [U 
[Till], and strongly consistent estimators (i.e. convergent almost surely), based on {dj^k){j,k)&i?i the 
discrete wavelet transform of Ifsm, have been proposed; notice that the latter estimators of H do not 
require that a to be known. Throughout our paper, for all (j. A;) G Z^, the wavelet coefficient dj^k is 
defined as, 

dj^k = 2^' / Xj^,^(t)V(2^i - k)dt- (1.3) 

Jm. 

moreover, we only impose to the analyzing wavelet ij^ a very weak assumption: ip is an arbitrary 
real-valued non-vanishing continuous function with a compact support in [0, 1] and it has 2 vanishing 
moments i.e. 

ij{s)ds = / sij{s)ds = 0. (1.4) 



It is worth noticing that we do not need that ^2^/'^ip{2^ ■ —k) : {j,k) € 7?^ be an orthonormal 
wavelet basis for L^(M). 

In view of the fact that the problem of the estimation of H is now well understood, from now on we 
assume the latter parameter to be known. Our goal is to construct, by using the wavelet coefficients 
idj,k)o<k<2i ^ a strongly consistent (i.e. almost surely convergent when j — )• +cxd) estimator aj of the 
parameter a. Let us outline the main ideas which lead to this estimator. 

• The starting point, is a result of [8j, according to which, with probability 1, the quantity 
H — 1/a, is the critical uniform Holder exponent of the sample paths of X}{,a over any arbitrary 
compact interval and in particular the interval [0, 1]; more precisely, one has, almost surely for 
all arbitrarily small rj > 0, 

j \XH,aiti) - XH,ait2,^ . 
sup < ; ,„ , > < OO (1.5) 

tut Jo,i]\ \ti-t2\"-y--^ ^ 

and 

tj:.m\ |ti-t2|^-V"+'^ 

• Next, let us set. 



■ fXH^aih) - XH,ait2)\ 1 

sup < — n — . } = (1-6) 



Dj = max \djk\- (1-7) 

0<k<23 



In view of the fact that the wavelet tp has a first vanishing moment, one can derive from (jl.5 
that, almost surely, for all arbitrarily small e > 0, 



limsup { 2^(^-1/"-^) I < OO. (1.8) 
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• Notice that, since we do not impose to "0 to be a continuously differentiable function and to 
|2J/V(2-'' • -k) : (j, fc) € Z2} to for m an orthonormal wavelet basis for L^(M), a priori it is not 
at all clear that (jl.6p . implies that, almost surely for all arbitrarily small e > 0, 



limsup |2-'(^-i/"+^)Z)A = oo. (1.9) 

Yet, by making use of some specific properties of Ifsm as well as the fact that -0 is compactly 
supported, we will be able to show that, a result stronger than ()1.9p holds; namely, one has 
almost surely, for all arbitrarily small e > 0, 

liminf (2^(-f^-i/°+^)i:>,- j = oo. (1.10) 

• Finally, combining (jl.Sp with (jl.lOp . one can get the following theorem, which is our main 
result. 

Theorem 1.1. For each j G N, one set, 

1 

aj 

where Dj is defined in |i. Tf ). Then, one has 

^ a.s. 

aj > a. 



log(£^ 
jlog(2)' 

almost surely. 



2 Proofs 

2.1 Proof of Relation ([HI]) 

The proof is standard in the wavelet setting, we give it for the sake of completeness. Let Cl be an 
event of probability 1 on which Relation (jl.Sp holds and let uj £ Cl be arbitrary and fixed. Assume 
that e > is arbitrary and fixed and denote by C{uj) the finite quantity defined as, 

C{u:):= sup <^ ^ ,H-i/a-e (■ (2.1) 

On the other hand, notice that ()1.3p . ()1.4p and the fact that, 

suppV' C [0,1], (2.2) 
imply that, for all (j, K) E Z+ x Z+ satisfying < < 2-', one has, 

r.(fc+l)2-J 



/ \xH,c.itM-XH,o.ik2-^M}m'i-'^)di. (2.3) 
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Next, combining (j2.3|) with (j2.1|) . one gets 



-(fc+l)2-J 



< ||^||io.(K)CM2W |t-A;2-^|^~^/"-^dt 

< ||^|U..(K)C(u;)2-^-(^-i/"-^), 



which proves that (jl.Sp is satisfied. □ 



2.2 Proof of Relation ( ITTOD 



Let us first recall that in [2j, a nice stochastic integral representation of the wavelet coefficients dj^k 
has been obtained, namely one has almost surely that 

d,, fc = 2-^(^-1/") [ ^hA'^'s - ids) , (2.4) 

JR 



^hA^)= / (y-x);-^/Xy)dy= / (y-x);-^/"V'(y)dy; (2.5) 



where ^H,a is the real-valued continuous function defined for each x € M, as, 

;+ ~^^°'tl^{y)dy = / 
Jo 

notice that the last equality results from (j2.2p . 

Proposition 2.1. T/ie function ^H,a satisfies the following two nice properties: 

(i) one has, 

supp$H,a ^ (-00,1]; (2.6) 

(ii) there is a constant ci > such for all x ^ [ — 00, 1], 

|^H,a(x)|<ci(l + N)-('+^/"-^^ (2.7) 

Proof of Proposition [KT[ Part (i) is a straightforward consequence of (|2.5|) and (|1.2|) . Let us show 
that Part (ii) holds. First observe that, (|2.5p easily implies that, 

sup |(l + |x|)'+'/""^|$j^,„(x)|| <4||V'||loo(k) <cx). (2.8) 

Let us now suppose that x < — 1. We denote by the primitive of ip, defined for all z € M, as 

i;(-^){z) = r i,{y)dy. 



Observe that ()2.2p and ()1.4p entail that the continuous function ij)^ has a compact support included 
in [0, 1]. We denote by -i/^^"^) the primitive of defined for all z € M, as 

V(-2)(z)= r ^l^^-^\y)dy. 
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Observe that supp';/''^ C [0, 1] and (|1.4|) entail that the continuous function -0^ has a compact 
support included in [0, 1]; therefore integrating two times by parts in ()2.5p . we obtain 

^hA^) = {H- l/a){H - l/a - 1) f\y - x)''-^l^-^^^-^\y)dy. (2.9) 

Jo 

Next, using ()2.9p and the inequalities: for all y G [0, 1], y — x > |x| > 2~^(l + \x\), it follows that, 

\^hA^)\ < 2^+^/"-^||^(-^)||loo(m)(1 + \x\f-'/''-\ (2.10) 
Finally, combining ()2.8p with ()2.10p . we get Part (ii) of the proposition. □ 
A straightforward consequence of ()2.4p and Part (i) of Proposition 12. H is that, 

d,,k = 2-^(^-1/") Z^'^'^' $H,a(2^s - A;)Z„ (0 ds). (2.11) 

^ — oo 



Let us now introduce some additional notations. We assume that 5 G (0, 1/3) is arbitrary and 
fixed. For all j G Z+, we define the positive integer ej as, 

ej := [2^'^], (2.12) 

where [•] is the integer part function. Then, for any integer / such that 

< / < [2J'(i-^)] - 1, (2.13) 

we set 

^(/e,+l)2-J- 

Gj,ie, := / ^hA'^'s - /e,)Z„ (() ds), (2.14) 

and 

((Z_l)ej.+l)2-i 

Rj,ie,:= $i^,„(2Js-/ej)Z„(()ds). (2.15) 

J —OD 

Thus, in view of (j2.1ip . the wavelet coefficient dj,zej can be expressed as. 

Now, our goal will be to derive the following two lemmas which respectively provide lower and 
upper asymptotic estimates for maxQ<;^[2J(i-'')] |G'j,Zej | and maxg^^^pjCi-*)] l-^j.^e^l- 

Lemma 2.1. One has, almost surely 

liminfi2^'^ max |G,;ell>l- (2.17) 

j^+oo [ 0<«<[2^(l-«)] ' J 

Lemma 2.2. One has, almost surely 



limsup<^2-'^ max \R.i^\\={). (2.18) 



5 



The proof of Lemma 1 2 . 1 1 mainly relies on the following two results. 

Lemma 2.3. (see e.g. l^) Let Y be an arbitrary symmetric a-stable random variable with a non- 
vanishing scale parameter ||5^||q, then for any real number t > ||^||a, one has, 

cslli^ll^r" < P(|y| >t)< C2||r||St-", (2.19) 

where C2 and C3 are two positive constants only depending on a. 

Lemma 2.4. For each fixed j € Z^-, {Gj^i^j '■ < I < [2-'(-'^^'^)] — 1} is a sequence of identically dis- 
tributed independent symmetric a-stable random variables whose scale parameters, denoted \\Gj^ie \\a, 
satisfy for all I, 

= 2-^' / l^HA^Tdx. (2.20) 

Proof of Lemma \2.4\ The independence of these symmetric a-stable random variables is a straight- 
forward consequence of the fact that they are defined (see ()2.14p ) through stable stochastic integrals 
over disjoint intervals. In order to show that they are identically distributed it is sufficient to prove 
that (|2.20p holds for each /. Using a standard property of stable stochastic integrals (see e.g. [S]) 
and ()2.14p . one gets 

\\Gj,ieX= / I'^nA'^'s-lejTds; 
then the change of variable u = 2^ s — lej allows to obtain ()2.20p . □ 



Now, we are in position to prove Lemma l2. 11 

Proof of Lemma \2.1l Let j € be arbitrary and fixed. Using the fact that {Gj^iej : < / < 
[2j'(i~'5)] — 1} is a sequence of independent identically distributed random variables (see Lemma . 
one gets, 

[2i(i-5)]-l 

Pf max \Gj^iA<2-^'^) = TT p(|G,- ■ | < 2"^'^ 



1=0 



P( |G,nl < 2-J'^'l =( 1 -PflG.nl > 2--''^n • (2.21) 



Observe that, in view of (|2.2Up in which one takes / = and in view of the assumption that 5 G (0, 1/3), 
there exist a positive constant C4 and a positive integer jo such that, one has, 

l/a 



2-^^>||G,ye,||a = 2-^/M / l^nA^Tdx] > ^'"2-^/°, (2.22) 



•^3 



for all integers j and I satisfying j > jo and < / < [2^^^ ^^]; notice that the last inequality in ()2.22p . 
follows from the fact that we have chosen jo, such that for every j > jo. 



[' \^HAxTdx>2-^ C \^HAxTdx, 

Jl—e-i J —00 
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and the last integral is positive since ^H,a is a non-vanishing function (this is a consequence of our 
assumptions on ip). Also notice that, one can suppose that C4 G (0,03^) (the positive constant C3 
has been introduced in Lemma l2.3p . Next, it follows from (j2.2ip . from the first inequality in ()2.19p 
in which one t = 2~^~^ , and from ()2.22p . that, for all integer j > jo, 

P( max |G,-,e,|<2-^-^)<(l-C52-^(i-2^))' (2.23) 

where the constant C5 := C3C4 G (0,1). Then, ()2.23p . the fact that b G (0,1/3), and standard 
computations, allow to show that, 

+00 

Vpf max |G,ie,| < 2-^'^') < 00; 

J=JO 

thus, applying the Borel-Cantelli Lemma, one gets ()2.17p . □ 
The proof of Lemma 12.21 mainly relies on the following result as well as on Lemma [ 



Lemma 2.5. For all non-negative integers j and I such that I < [2^^], the scale parameter \\Rj^ie \\a 
of the symmetric a-stable random variable Rj^ie^ (see 112.15]) ) satisfies, 

\\Rj,ie^\a = 2"'' /' l^nA^Tdx < ca-M2S+i/a-5H)^ (2.24) 

where cg is a positive constant non depending on j and I. 

Proof of Lemma \2.5l The equality in ()2.24p can be obtained by using ()2.15p and the arguments which 
have allowed to derive (j2.20p . Let us show that the inequality in (|2.24p holds; there is no restriction 
to assume that j > 6~^. Using ()2.7p and ()2.12p . one has, 

2-^ f |$H,a(x)r(ix < 05^2-^' f ''(1 

J —oo J —oo 
< C?2-^- / (1 + x)-'°-'+"^dx = C? ^\ \r. < C62-^°(25+l/«-5^^)^ 

J2JS-2 a{2- H) 



where the constant 

2a{2-H) 



C6 := c° 



^a{2-Hy 

□ 



Now, we are in position to prove Lemma l2.2i 



Proof of Lemma \2.2[ First, observe that in view of the assumption that 5 G (0, 1/3), one has for a 
fixed arbitrarily small rj > 0, 

+ ^ < 2<5 + 1/a - 6H; 



a 
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therefore, it follows from Lemma l2.5t that there exists a positive integer j'l, such that for all integers 
j and I, satisfying j > ji and < / < [2^^-^~^^, one has, 

Thus, we are allowed to apply the second inequality in (j2.19|) . in the case where Y = Rj^ie^ and 

■ ( 2S + r, \ 

t = 2 ■'y °' /.As a consequence, we obtain that, for all j > ji, 

Pf max \Ru^.\>2-^i^))< V pf . I > 2"^(^) ) 
Vo<«<[2^(i-*)]' ^ ~ ^ ^ ^ 

[2i(i-5)]_l 

<c^2^i^s+n) J2 ||i?,,,e,||° <C72-^'"(25+i/"-^^^)+i(i+^+'?), (2.25) 

1=0 

where the last inequality results from (|2.24p and the constant cy := C2Cq. Assume that 6{a — 1) > ij, 
then one has, 

a{26 + 1/a - 6H) > a{6 + 1/a) = a5 + 1 > 1 + 6 + i]. 
Therefore, it follows from ()2.25p that, 

+ 00 

Vpf max |i?,ze I > 2"^A^)) < oo; 
^ Vo<Z<[2J(i-«)l ' J 
3=31 - ^ ^ 

thus, applying the Borel-Cantelli Lemma, one gets (j2.18p . □ 

Remark 2.1. Our proof s of Lemmas \2. 1\ and \2.^ only allow to derive that Relations IjS.lT^ and i2.18\) 

hold on some event of probability 1, denoted by Qs, since it a priori depends on 6 ^ (0, 1/3). Yet, one 
can easily show that these two relations also hold, for every real number 5 G (0, 1/3), on an event of 
probability 1 which does not depend on 5, namely the event n5GQn(o i/3) 

Now, we are in position to prove Relation (jl.lOp . 
Assume that e is a fixed arbitrarily small positive real number and that 5 G (0, 1/3) is such that, 

e/2 = 25/a. (2.26) 

Next observe that (j2.26p , (|1.7p , (|2.16p and the triangle inequality, imply that for all j G Z+ , 

2i(H-i/a+e/2)£) > 2Jlr max |G. i^. + R, leA 

0<i<[2J(i-*)] ' 

> 2-' Q max \Gjie\—2-'a max LR,/e ; 

0<i<[2J(i-'5)] ^ 0<Z<[2i(i-5)] ^ 

therefore, one has that, 

liminf |2^(^-^/"+^/2)^ A (2.27) 

f 1 — 1 1 f 1 — 

> lim inf < 2-' ^ max \Gjie \ \ — lim sup < 2-' « max \Rj le ■ 

j^+oo [ 0<«<[2J(i-5)] ' J j^+oo I 0<R[2J(i-*)] ' 

Finally putting together, (l2T7l) . (I2T8]) . and ([2:26]) . one gets (fLTOD . 
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2.3 Proof of Theorem [TTTl 



Relations (jl.Sp and (jl.lOp imply that there is an event of probability 1 such that each a; € 17* 
satisfies the following property: for all arbitrarily small e > 0, there are two finite positive constants 
A = A{io,e) and B = B{uj,e), and there exists j2 = j2{uj,e) G Z+, such that, one has for all integer 
J > J2, 

This entails that, 

-H + l/a-e<lim inf | l^liM^ \ < Ih, sup | i^liM^ 1 < + l/« + 
i^+oo 1^ jlog(2) J j^+oo I jlog(2) J 

Then letting e goes to zero, one gets that, 

jlog(2) j ^ 

□ 
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